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ABSTRACT 

The  turbulent  mass  flux,  or  equivalently  the  fluctuating  Favre  velocity  mean,  appears  in 
the  first  and  second  moment  equations  of  compressible  k  —  e  and  Reynolds  stress  closures. 
Mathematically  it  is  the  difference  between  the  unweighted  and  density- weighted  averages  of 
the  velocity  field  and  is  therefore  a  measure  of  the  effects  of  compressibility  through  variations 
in  density.  It  appears  to  be  fundamental  to  an  inhomogeneous  compressible  turbulence, 
in  which  it  characterizes  the  effects  of  the  mean  density  gradients,  in  the  same  way  the 
anisotropy^  tensor  characterizes  the  effects  of  the  mean  velocity  gradients.  An  evolution 
equation  for  the  turbulent  mass  flux  is  derived.  A  truncation  of  this  equation  produces  an 
algebraic  expression  for  the  mass  flux.  The  mass  flux  is  found  to  be  proportional  to  the  mean 
density  gradients  with  a  tensor  eddy-viscosity  that  depends  on  both  the  mean  d^ormation 
and  the  Reynolds  stresses.  The  model  is  tested  in  a  wall  bounded  DNS  at  Mach  4.5  with 
notable  results. 
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tract  No.  NASl-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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1.  Introduction 


This  article  presents  a  derivation  of  a  representation  for  the  non-zero  first-moment  of 
the  iBuctuating  velocity  field,  the  time  average  of  the  fluctuating  component  of  the  Favre 
velocity,  <  v,  >.  Mathematically  <  v,  >  represents  the  difference  between  unweighted  and 
density  weighted  averages  of  the  velocity  field  and  is  therefore  a  measure  of  the  effects  of 
compressibility  through  variations  in  density.  It  plays  an  imp>ortant  role  in  parameterizing 
the  anisotropic  effects  of  compressibility  associated  with  the  mean  dilatation  and  gradients  in 
the  mean  velocity  and  density.  Experimentally  it  is  an  important  quantity  that  allows  Favre- 
averaged  numerical  results  to  be  related  to  time-averaged  experimental  results.  The  need  to 
consider  this  quantity  is  motivated  by  its  frequent  contributions  to  the  first  and  second-order 
moment  equations  in  two-equation  k—e  type  turbulence  closures  as  weU  as  in  Reynolds  stress 
closures.  In  the  mean  momentum  equations  the  mass  flux  makes  a  contribution  to  the  viscous 
terms.  In  the  mean  energy  equations  the  mass  flux  makes  a  contribution  to  the  viscous,  the 
pressure  work,  and  the  pressure  flux  terms.  In  the  Reynolds  stress  equations  the  viscous 
terms  appear  naturally  in  Reynolds  variables  while  the  problem  is  posed  in  Favre  variables. 
In  the  process  of  splitting  the  viscous  terms  into  the  viscous  transport  terms,  earned  in  Favre 
variables,  and  the  dissipation  terms,  carried  in  Reynolds  variables,  important  contributions 
from  the  mass  flux  appear.  The  accurate  accounting  of  these  terms  is  important  for  any 
consistent  near  wall  modeling  and  the  retention  of  the  mass  flux  terms  is  important  in 
complex  compressible  turbulent  flows.  These  contributions  have  been  investigated  in  detail 
in  Ristorcelli  (1993).  The  mass  flux  also  determines  the  importance  of  two  production 
mechanisms  one  due  to  the  acceleration  of  the  mean  flow  and  the  other  due  to  viscous  effects 
associated  with  the  Favre  fluctuation  mean. 

Many  of  these  contributions  are  neglected  in  turbulence  closure  models.  This  is  a  result  of 
assuming  that  the  Favre  mean  velocities  are  suitable  approximations  to  the  Reynolds  mean 
velocities.  This  approximation  is  not  appropriate  in  complex  flows  of  practical  interest.  The 
retention  of  the  mass  flux  terms  will  be  necessary  in  complex  compressible  turbulent  flows: 
these  include  flows  in  which  there  are  mean  density  gradients  due  to  large  Mach  number, 
combustion,  separation  or  reattachment  (inflection  points),  cold  wall  boundary  conditions, 
mean  dilatation,  shocks,  adverse  pressure  gradients,  or  strong  streamwise  accelerations.  Even 


in  this  nominally  simple  compressible  flow,  such  as  a  supersonic  wall  bounded  boundary  layer 
which  has  a  four-fold  variation  of  the  mean  density  over  the  width  of  the  boundary  layer,  the 
mass  flux  is  not  negligible.  Dinavahi  et  al.  (1993b),  in  a  Mach  4.5  wall  botmded  DNS,  has  n 


found  that  the  cross-stream  Favre  mean  and  Reynolds  mean  velocities  have  different  signs 


attesting  to  the  fact  that  the  mass  flux  is  not  small  with  respect  to  the  mean  velocities. 


Others  have  recognized  the  importance  of  the  man  flux  and  several  models  have  been 
proposed.  Taulbee  and  VanOsdol  (1991)  have  derived  a  modeled  equation  for  the  man 
flux.  In  thdr  equation  they  keep  the  correlations  with  the  surface  forces  which  are  modeled 
nsuming  a  homogeneous  turbulence  and  the  validity  of  Morkovin’s  hypothesis.  In  the  present 
asymptotic  derivation  these  terms  scale  with  the  density  intensity  and  are  found  to  be  of 
higher  order;  the  present  asymptotic  derivation  of  the  transport  equation  for  the  man  flux 
keeps  only  the  zeroeth  terms  to  keep  the  model  simple  and  to  avoid  the  Ion  of  accuracy 
associated  with  the  individual  sq>proximati<ms  made  in  many  models.  They  also  use  a 
gradient  transfer  anumption  for  the  turbulent  diffusion.  Due  to  the  different  manipulation 
to  obtain  an  equation  for  the  the  Favre  fluctuation  equation,  the  turbulent  diffusion  is  fotmd 
to  be  proportional  to  the*  difference  between  the  Favre  and  Reynolds  averaged  Reynolds 
stresses;  as  this  is  a  higher  order  term,  scaling  as  the  density  intensity,  there  is  no  need  to 
model  it.  The  Taulbee  and  VanOsdol  model  requires  the  solution  of  two  modeled  differential 
equations,  one  for  the  mass  flux  and  one  for  the  density  variance  which  «q>pears  as  a  source 
term  in  their  modeled  mass  flux  equation.  In  the  present  derivation  of  the  evolutimi  equation 
for  the  mass  flux  only  one  unknown,  the  covariance  with  the  fluctuating  dilatation,  requires 
modeling.  There  is  no  need  for  a  separate  equation  for  the  density  variance. 

Zeman  and  Coleman  (1991)  have  also  proposed  a  mass  flux  model.  Their  modeled  equa¬ 
tion,  which  has  been  tested  in  the  turbulence  through  a  shod(  simulations  of  Lee  (1992)  is 
very  similar  to  the  one  derived  in  this  article.  They  also  propose  an  algebraic  expression  for 
the  mass  flux  to  which  the  present  model  simplifies  to,  in  the  limit  of  negligible  mean  defor¬ 
mation.  Otir  work  has  shown  that  the  inclusion  of  the  mean  velocity  gradients  is  essential 
to  capturing  the  near  wall  maxima  of  the  mass  flux.  After  all,  the  mean  velocity  gradients 
are  a  major  portion  of  the  production  terms  of  the  mass  flux  in  its  transport  equation. 

Rubesin  (1990)  has  also  proposed  a  mass  flux  model.  It  assumes  that  1)  the  fluctuations 
obey  a  polytropic  gas  law  2)  the  specific  heats  are  constant  allowing  the  fluctuating  density 
to  be  written  in  terms  of  the  fluctuating  enthalpy  and  that  3)  the  fluctuating  enthdpy  can 
be  related  to  the  mean  enthalpy  using  a  gradient  transfer  hypothesis.  The  Rubesin  model 
requires  the  poly  tropic  index  as  an  input.  Dinavahi  et  al.  (1993)  and  Ristorcelli  et  at.  (1993) 
in  a  temporal  DNS  have  shown  that  the  polytropic  index  varies  substantially  over  the  width 
of  the  turbulent  boundary  layer.  The  Rubesin  model  also  predicts  a  mass  flux  only  when 
there  is  a  heat  flux,  while  the  present  model  derived  from  the  exact  evolution  equation, 
predicts  a  mass  flux  whenever  there  are  mean  density  gradients. 

The  present  model  for  the  mass  flux  starts  with  the  exact  evolution  equation  for  the  mass 
flux.  An  equation  for  the  fluctuating  Favre  mean  is  then  developed  in  a  power  series  in  the 
fluctuating  density  intensity.  To  zeroeth  order  there  is  only  one  unknown  correlation  and 
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no  unknown  correlations  with  viscous  and  pressure  terms  appear.  The  evolution  equation  is 
simple  enough  to  carry  as  an  additional  differential  equation  in  turbulence  simulations,  as 
Zeman  and  Cole  (1991)  have  proposed.  Nonetheless  an  algebraic  truncation  of  the  equation 
is  derived  as  a  further  simplification  of  the  problem  applicable  to  most  compressible  flows  of 
engineering  interest.  The  algebraic  tnmcation,  similar  to  that  used  in  algebraic  stress  models, 
assumes  a  structural  equilibrium  which  relates  the  material  derivative  in  the  fluctuating 
Favre  mean  equation  to  the  production  and  dissipation  in  the  kinetic  energy  equation.  The 
truncation  produces  a  set  of  three  coupled  algebruc  equations,  of  the  form  Aij  <  Vj  >=  6,. 
Application  of  the  Cayley- Hamilton  theorem  produces  an  explicit  closed  form  expression 
for  the  <  Vi  >.  The  <  v,  >  are  found  to  be  proportional  to  the  density  gradients  with 
an  eddy-viscosity  tensor  dependent  on  the  Reynolds  stress  and  the  mean  deformation.  The 
fluctuating  Favre  mean  is  then  related  back  to  mass  flux  using  the  well  known  relation 
between  the  two  quantities. 

This  article  is  organized  in  the  following  manner.  After  motivating  the  investigation  in 
section  two,  section  three  describes  the  derii^tion  of  an  evolution  equation  for  the  mass 
flux.  In  the  following  section  an  algebraic  model  for  the  mass  flux  is  obtained.  The  general 
expression  for  the  <  >  is  then  specialized  to  several  simple  mean  flows  in  order  to  highlight 

the  physics.  It  is  found  that,  in  the  limit  of  isotropic  turbulence  with  nqgli^ble  mean  velocity 
gradients,  the  derived  expression  reduces  to  the  usual  scalar  eddy- viscosity  form  derived  using 
a  gradient  transfer  assmnption.  The  model  is  tested  in  the  Ma  »  4.5  waU  bounded  DNS  of 
Dinavahi  and  Pruett  (1992). 

2.  Preliminary  exposition 

In  general,  upper  case  letters  will  be  used  to  denote  mean  quantities  except  in  the  case  of 
the  mean  density,  <  p  >,  since  p  which  has  no  convenient  upper  case  form.  The  averaging 
operation  is  indicated  using  the  angle  brackets  for  time  means,  <  ViVj  >,  and  the  curly 
brackets  for  the  density-weighted  or  Favre  mean,  where  <  p  >  =<  p^Vi^j  > 

and  the  asterisk  denotes  the  full  field,  p*  =<  p  >  +  pf.  The  dependent  variables  are 
decomposed  according  to 

n*  =  Ui  +  Ui  where  <  u<  >=  0 

ti*  =  -1-  where  {v<}  =  0 

p*  =  <  p  >  +  p*  where  <  p'  >=  0 

p*  =  P  +  p  where  <  p  >=  0 

T*  =  T  +  0  where  =  0 

As  both  the  Reynolds  and  the  Favre  velocities  appear  naturally  in  the  evolution  equations 
for  a  compressible  turbulence  it  is  necessary  to  carry  both  the  Favre  and  the  Reynolds 
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decc»Dipo«ti<Hi8  ol  the  velocity  field.  They  are  related  by 


Ui=  K  +  <  > 

m  =  Vi  —  <  Vi  >  . 

The  fluctuating  Favre  mean  quantifies  the  difference  between  the  Favre-mean  and  the  Reynolds- 
meui  velocities,  Vi  and  l/j,  as  well  as  the  difference  between  the  instantaneous  fluctuating 
portions  of  these  two  fields.  Note  that  because  of  the  definition  of  the  Favre-average  of  the 
Favre-deviation,  <  p*Vi  >=<  p  >  {«,}  =  0, 

—  <  P  ><  Vi  >=<  pVi  >, 


the  time  average  of  the  fluctuating  Favre  velocity  and  the  mass  flux  are  equivalent  quantities 
(apart  from  a  scaling  by  the  local  mean  density).  Because  of  the  peculiarities  of  the  density- 
weighted  averaging  operation,  a  second-order  statistic,  <  pVi  >,  can  be  expressed  as  the 
product  of  two  first-order  statistics,  —  <  p  ><  u,-  >.  The  two  phrases  mass  flux  and  Favre 
fluctuation  mean  will  be  used  interchangeably.  The  primes  on  the  fluctuating  density  have 
been  dropped. 

As  Ui  =  Vi  +  <  Vi  >,  the  <  Vi  >  quantify  the  difference  between  the  unweighted  or 
Reynolds  mean,  and  the  density- weighted  mean,  V{,  and  represent  the  effects  of  compress¬ 
ibility  through  variations  in  density.  Data  from  Ma  =  4.5  DNS  computations  of  Dinavahi 
and  Pruett  (1993)  in  unidirectional  developing  wall  bounded  flow  indicate  that  the  approx¬ 
imation  of  Ui  a  Vi  in  the  wall  region  is  inadequate.  In  this  flow,  in  which  Mt  ^  0.3  and 
there  is  a  four-fold  variation  of  the  mean  doisity  over  the  boundary  layer.  In  data  taken 
from  that  simulation,  shown  in  Figure  1,  it  was  unexpectedly  found  that  <  V2  >  »  larger 
than  either  and  V3.  It  is  large  enough  to  cause  U2  and  V2  to  have  different  signs.  This 
is  an  indication  that  the  net  fluid  particle  transport  and  the  net  momentum  transport  are 
in  opposite  directions.  The  point  is  that  this  is  a  nominally  simple  flow,  in  comparison  to 
those  of  practical  interest,  no  inflection  points,  no  change  of  geometry,  no  substantial  heat 
transfer,  no  cold  wall  boundary  conditions  with  the  concomitant  change  in  sign  of  the  mean 
density  gradient,  in  which  the  approximation  U2  —  V2  was  expected  to  be  adequate  and  they 
are  not  even  of  the  same  sign. 

In  compwing  experimental  data  and  computational  results  the  mass  flux  plays  a  role  in 
relating  the  Reynolds  stresses  in  Favre,  Uj,  variables  and  Reynolds,  u,-,  variables: 

<  ViVj  >=  <  UiUj  >  +  <Vi  ><  Vj  >  . 

The  moments  involving  u,-  are  experimentally  measured  while  those  involving  v,-  come  from 
the  calculations.  As  it  is  a  vector  it  describes  the  anisotropic  effects  compressibility  has 
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on  the  turbulence.  An  anisotropy  tensor  based  on  the  Reynolds  variables  is  defined  as 
bij  =<  UiUj  >  /  <  UpUp  >  —  IfZSij.  A  similar  anisotropy  tensor  using  the  Favre  variables 
can  also  be  defined:  =<  ViVj  >  /  <  VpVp  >  —  An  energy  weighted  deviation  of  the 

anisotropy  tensor  from  its  density-weighted  equivalent  is  given  solely  in  terms  of  the  Favre 
fluctuation  mean: 


<  VpVp  >  b(j—  <  tipUp  >  bfj  =<  Vi  ><  Vj  >  —1/3  <  Vg  ><  Vg  >  Sij, 


Note  that  there  are  only  three  independent  quantities  <  v«  >.  ^FVom  a  heuristic  point  of 
view  this  is  pleasantly  consistent:  the  effect  of  mean  flow  gradients,  K',ii  is  parameterized  by 
the  six  components  of  the  anisotropy  tensor  while  the  effect  of  the  mean  density  gradients, 
the  vector  <  p  >,i,  is  parameterized  by  the  three  components  of  the  mass  flux. 

There  are  some  interesting  properties  of  the  mass  flux  that  can  be  surmised  from  the 
above  relationships.  The  most  striking,  and  this  is  a  rigorous  result,  is  that  in  1)  ui  isotropic 
turbulence  and  or  in  2)  a  statistically  stationary  homogeneous  turbulence  with  mean  velocity 
gradients  and  with  no  mean  density  gradients,  <  v,-  >=  0  and  the  Reynolds  and  Favre 
variables  are  equivalent: 

Vi 


u,  = 
<  ViVj  >— 


Vi 

<  UiUj  > 


Similar  results  hold  for  relationships  between  the  various  moments  of  Ui  and  Vi  as  can  be 
easily  derived.  These  results  come  from  the  following  two  facts:  1)  in  an  isotropic  field  all 
vector  statistics  are  zero;  2)  in  a  statistically  stationary  homogeneous  field,  whose  directional 
characteristics  are  solely  determined  by  the  mean  velocity  gradients,  which  are  invariant  to 
coordinate  reflection,  all  quantities  not  invariant  to  coordinate  reflection  are  zero.  This 
has  been  recognized,  in  the  context  of  compressible  turbulence,  by  Blaisdell  (1991).  In 
short,  in  isotropic  or  homogeneous  turbulence  without  mean  density  gradients,  there  is  no 
difference  between  the  problem  posed  in  Reynolds  or  Favre  variables.  This  is  an  important 
and  serious  issue  affecting  the  vididity  of  conclusions  about  the  performance  of  compressible 
turbulence  models  which  have  been  developed  and  tested  in  homogeneous  or  isotropic  flows. 
On  the  other  hand  it  suggests  the  appropriateness  of  the  incompressible  turbulence  modeling 
framework  in  building  models  for  the  compressible  flow  as  they  are  consistent  in  the  isotropic 
and  homogeneous  limit,  for  arbitrary  turbulent  Mach  number. 
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The  second  ooKmiesit  equaticms  for  a  compmsibleflow  are  written,  without  approximation 
and  alter  some  manipulation,  as 

DfDt  (<  />  >  {».«,})  =  -<p>  {w.Vp}  -<p>  +  Ily  +2/3  <  jw*,*  >  Sij 

-[<  >  Spj+  <  pvj‘>  6ip+  <p>  {viv^v,}-  <  >  -  <  Vi<r)^  >],p 

+  <  Wj  >  (-P,i  +  E<*,*  +  <  <T<*  >,*  1+  <  Vi  >  [-P„  +  E,fc,*  +  <  <Tik  >yk 

-  <  u,,parj  >  -  <Ui,pff^> 

where  the  mean  momentiun equations  have  been  used  and  (Tij  =<  p  >  [vi,j  +Vj,i  — 2/3v,„ 

Ei,  =<  p>  {l^,j+Vj,<-2/3V,„^ij]  and  p  >  [tti,j +Uj„ -2/3«„,^<j].  The  form  of 

the  equations  above  reflects  tbe  following  manipulations:  1)  The  deviatoric  part  of  the 
pressure-stridn  correlation  is  defined  as  n«i  ===<  P(Vui  +Vit« )  >  —2/3  <  pvfc,fc  >  6^  and  2) 
the  identity  Vj  =  u,-  +  <  v,  >  has  been  used  to  rewrite  the  transport  terms  in  v,  vui- 
ables  while  keeping  the  dissipation  terms  in  u,-  variables.  In  the  equations  for  the  Favre- 
averaged  Reynolds  stress  the  terms  arising  from  surface  forces  appear  naturally  in  (f/«,Ui) 
variables  while  the  problem  is  posed  in  (K',Vt)  variables.  In  recasting  the  Reynolds  variables 
terms  in  Favre  variables  the  mass  flux,  <  v«  >,  makes  several  different  contributions  to 
the  Reynolds  stress  equations  and,  of  course,  to  the  k  =  l/2{vjVj}  equation.  It  multiplies 
the  mean  flow  acceleration  which  is  a  new  turbulence  production  mechanism  important  in 
flows  with  strong  mean  pressure  gradients,  shocks  and  expansion  fans,  and  in  any  flows 
that  have  strong  streamwise  accelerations.  The  mass  flux  also  contributes  to  the  viscous 
diffusion  of  the  Reynolds  stresses  a  term  that  is  important  in  the  near  wall  region  which  is 
also  where  the  mass  flux  terms  are  important.  Note  that  =  crij—  <  (Xij  >  aUows  the 
viscous  transport  terms  to  be  recast  in  the  Favre  variables  and  that  mass  flux  terms  and 
their  derivatives  will  appear.  The  mass  flux  also  contributes  to  the  Reynolds  stress  equations 
through  the  pressure  flux  to  which  it  is  coupled  by  the  equation  of  state:  for  an  ideal  gas 
<  pvj  >=  P[<  pvj  ><  p  >~*  +  T"*].  In  the  adiabatic  case  the  pressure  flux  can 

be  written,  to  first  order,  as  <  pvj  >—  P'^  <  pVj  ><  p  >“*=  <?  <  pVj  >.  Results  from 
some  numerical  simulations  have  shown  that  the  pressure  and  density  fluctuations  of  the 
turbulence  passing  through  a  weak  shock  can  be  related  through  such  a  rule,  Lee  (1992). 

This  is  not  found  to  be  true  for  the  wall  bounded  flow  of  Dinavahi  and  Pruett  (1993)  as 
shown  in  Dinavahi  tt  al.  (1993).  Lee  has  also  found  that  the  pressure  flux  (as  well  as 
the  pressure-dilatation)  is  primarily  responsible  for  the  rapid  evolution  of  turbulent  kinetic 
energy  downstream  of  a  weak  shock. 

In  the  mean  momentum  and  mean  energy  equations  the  viscous  terms  appear  in  Reynolds 
variables,  Ey(l/)  =<  p  >  +f/;M  — 2/3(/,,,tfij].  When  the  problem  is  recast  in  Favre 

variables,  the  viscous  terms  become  functions  of  the  Favre  mean  velocity  and  the  Favre 
fluctuation  mean.  It  is  typical  to  approximate  (/,  to  close  the  equation.  This  involves 
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neglecting  <  u,  >  as  I/,  =  K  4-  <  u,  >.  It  is  clear  that  this  approximation  is  only  valid 
when  <  Vi  >  and  its  gradients  are  negligible.  The  data  of  Dinavahi  et  al.  (1993b)  indicates 
that  this  approximation  is  a  poor  one  in  the  wall  bounded  flow.  In  fact,  in  some  portions 
of  the  turbulent  boundary  layer,  the  mass  fluxes’  contribution  to  the  viscous  terms  is  as 
high  as  25%.  Figure  2,  taken  from  Ristorcelli  tt  al.  (1993),  shows  the  second  cross-stresun 
derivatives  of  the  f/,-,  Vi  and  <  Wi  >.  The  mass  flux  terms  also  contribute  to  the  pressure 
and  viscous  work  terms  in  the  mean  energy  equation. 

It  is  clear,  given  the  number  of  times  it  occurs  in  the  moment  evolution  equations,  that  an 
expression  for  the  mass  flux  for  general  compressible  turbulent  flows  of  aerodynamic  interest, 
is  necessary.  An  evolution  equation  and  a  model  for  <  v,-  >  are  important:  1)  to  be  able  to 
estimate  the  importance  of  <  Vj  >  in  different  flows,  2)  to  know  what  to  do  about  it  when 
it  is  important,  and  3)  to  be  able  to  relate  experimental  values  to  computationsd  results. 

3.  An  evolution  equation  for  the  Favre  velocity  perturbation 

Consider  the  evolution  equations  for  the  total  velocity  and  density  fields: 

p\i  +  =  0 

(/»*<)><  +  (/“JOv  =  -p;  +  pV*  +  <^iij 

where  (T^j  =  +u^,i  -2/3u;„  To  a  very  good  approximation  the  viscosity  is  in¬ 

dependent  of  density:  it  will  be  taken  to  be  equal  to  its  local  mean  value  and  correlations 
between  the  viscosity  and  velocity  will  be  considered  as  higher  order  effects  and  neglected. 
The  evolution  equation  for  the  fluctuations  around  the  Favre-mean  momentum  are  obtained 
by  subtracting  the  evolution  equation  for  the  mean  momentum  <  p  >  Vi  from  the  equation 
for  p*u*  to  obtain  an  equation  for  p*Vi.  As  <  p*v,  >=  0,  a  straightforward  time-average  of 
this  equation  does  not  produce  any  results.  The  equation  is  rewritten  in  its  nonconservative 
instantaneous  form: 

P*Vi,t  +  p*(Vp  +  Vp)vi,p  +  {p*VpVi),p  -  <  p^ViVp  >,p  +  2€ikpUkP*Vp  ^ 

-Pi  +  P*Si  +  -  plvjVi  Vi(p*Up),p 

where  =<  p  >  [u,-,j +u_,„ -2/3u„,  5ij]  and  =  Vi,f¥VpVi,p+2eikp^kVp  -  Fi.  The 
term  pLyjVi  reflects  the  coupling  between  the  fluctuating  density  and  the  mean  flow.  Divid¬ 
ing  by  the  total  density,  p*,  expanding  using  the  binomial  theorem,  and  averaging  pro¬ 
duces  an  evolution  equation  for  <  >  in  which  successive  terms  scale  as  where 

<T  =<  p^  >  /  <  p  is  the  normalized  density  variance. 

Keeping  only  lowest  order  terms  produces 

<  Wi  >,t  +  <Vp>=  -  <Vp>  Vi,p  -t-  {u.Vp}  <  P>,p<  p  >~^  +  <  ViVp,p  > 

+  ({t'pwj-  <  VpVi  >],p  -  <  fi>  -i-O  (y/ff) 
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an  evolution  equation  independent  of  complicating  correlations  with  the  pressure  and  viscous 
surface  forces.  Note  that  the  fact  that  <  ViVj  >=<  pViVj  ><  p  >“*  and  that  < 

p*ViVp  >=<  p  >  {ujUp}  has  been  used.  The  inhomogeneous  diffusion  term  [{v,Vj}—  <  ViVj  >] 
is  an  O  (>/<7)  term,  as  can  be  seen  by  the  data  presented  in  Dinavahi  et  al.  (1993),  and  can 
in  general  be  neglected.  In  a  homogeneous  turbulrace  it  is  zero,  of  course. 

This  very  simple  equation  for  <  Vi  >  results  from  the  fact  that,  in  the  Favre  setting, 
surface  forces  are  carried  using  the  Reynolds  decomposition  while  volume  forces  appear 
naturally  in  the  Favre  variables.  The  first-moment  of  the  fluctuating  surface  forces  (pressure, 
viscosity)  appearing  in  the  equation  for  <  Uj  >  are  zero  and  no  complicating  models  for  these 
terms  are  required.  This  combined  with  the  peculiarity  of  the  fluctuating  Favre  mean  allows 
<  pvi  >=  —  <  p  ><  Vi  >  and  leads  one  to  work  with  the  first-moment  form  <  v,  >  of  the 
second-moment  <  pvi>.  Thus  a  simple  evolution  equation  for  the  mass  flux  that  highlights 
the  zeroeth  order  effects  associated  with  the  volumetric  compressibility  while  relegating  the 
higher  order  effects  of  the  surface  forces  to  a  higher  order  equation  in  the  expansion  is 
obtained. 


4.  An  algebraic  expression  for  the  Favre-velocity  perturbation 

To  obtsun  the  mass  flux,  <  pVi  >==  —  <  p  X  Vi  >,  an  equation  for  the  Favre  fluctuation 
mean,  <  >,  with  only  one  unknown  term,  the  correlation  with  the  fluctuating  divergence 

<  ViVptp  >,  has  been  derived.  The  evolution  equation  obtained  for  the  Favre  fluctuation 
mean  is  simple  enough  to  carry  in  turbulence  simulations.  However  it  is  stiU  simpler  and 
less  expensive  to  carry  an  algebraic  expression.  This  is  now  derived. 

A  direct  algebraic  tnmcation  of  the  evolution  equation  will  describe  the  fixed  points  of  the 

<  Vi  >.  An  algebraic  truncation  following  the  procedure  used  in  algebraic  stress  models  will 

give  the  fixed  points  of  <  Vj  >  /{vpVpY^^,  This  is  done  by  assuming  a  structural  equilibrium 
of  the  form  DfDt  [<  u,-  >  =  0  aUowing  the  convective  derivatives,  DIDt  <Vi> 

to  be  expressed  in  terms  of  the  right  hand  side  of  the  evolution  equation  for  the  turbulence 


energy: 


D  <Vi>  D  f  ^  \  It 


which  allows  the  evolution  of  the  <  v,-  >  to  reflect  the  changes  in  the  energy  of  the  local 
turbulence  field.  Here  V,  e,  are  the  production  and  the  dissipation  in  the  turbulent  kinetic 
energy  equation  where  k  =  1  /2{vpVp}  is  the  specific  kinetic  energy. 

In  the  near  wall  re^on,  where  the  mass  flux  is  expected  to  be  the  most  important,  the 
flow  will  attain  a  structural  equilibrium  rapidly  and  such  an  approximation  will  be  adequate. 
Note  that  the  equality  V  =  e  corresponds  to  the  fixed  point  DfDt  <  v,  >=  0.  The  algebraic 
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fonn  of  the  evolution  equation  for  Favre  fluctuation  mean  is  now: 

<  Vi  >{V-  e)/k  =  -  <Vp>  Vi,p  <p>,p<  p  >■*  +  <  ViVp,p  >  . 

The  body  force  terms  and  the  Coriolis  terms  have  not  been  carried,  however  the  analysis 
can  be  carried  quite  easily  with  them  as  they  do  not  constitute  unknown  terms  that  require 
closure. 

It  remains  to  close  the  last  term  on  the  right  hand  side.  It  is  possible,  in  situations 
with  large  density  and  velocity  gradients,  to  neglect  the  correlation  with  the  fluctuating 
divergence.  This  is  equivalent  to  the  assumption  that  the  mean  flow  gradients  of  density  and 
velocity  are  large  and  set  the  balance  to  lowest  order.  It  can,  however,  be  shown  that  the 
correlation  with  the  fluctuating  divergence  scales  with  mean  flow  gradients  and  is  therefore 
not  negligible  in  a  general  flow.  Moreover  there  are  times  when  the  difference  between  the 
mean  production  terms  is  small  which  means  that  the  contribution  from  <  v,Vp,p  >  will  be 
important. 

The  correlation  with  the  divergence  will  be  represented  by  a  lineu'  relaxation  model. 
This  linear  relaxation  model  is  chosen  on  the  grounds  that  <  ViVp,p  >  and  <  v,  >  have  the 
same  tensorial  properties.  Both  belong  to  the  same  symmetry  groups,  satisfying  the  same 
reflectional  and  rotational  properties,  vanishing  in  isotropic  turbulence  and  in  an  equilibrium 
homogeneous  turbulence.  From  a  computational  point  of  view,  a  linear  relaxation  form  is 
desirable  as  it  avoids  the  possibility  of  a  singularity  in  the  inversion  of  the  velocity  gradient 
during  a  computation  and  is  consistent  with  realizability.  A  linear  relaxation  with  time  scale 
Td 

<  ViVp,p  >=  -  <  Vi  >  /Td 

is  chosen.  Zeman  and  Coleman  (1991)  have  sdso  used  a  linear  relaxation  with  acoustic  time 
scale  for  this  correlation.  The  time  scale,  Td,  in  the  model  for  <  ViVp,p  >=  —  <Vi>  /Td  may 
be  thought  of  as  a  dilatational  time  scale.  The  time  scale,  Td,  is  the  only  phenomenological 
parameter  assiuned  to  obtained  the  present  model  for  the  mass  flux.  Computations  using 
the  acoustic  time  scale,  Td  =  Mtk/e,  to  represent  the  dilatational  time  scale  have  been 
successful.  The  present  model  will  use  this  approximation  for  the  dilatational  time  scale. 
There  are  other  possibilities  though  at  this  time,  given  the  success  of  the  present  model, 
there  is  little  motivation  for  further  investigation. 

Substituting  for  the  unknown  correlation  with  the  fluctuating  divergence  in  the  algebraic 
truncation  of  the  evolution  equation  for  <  Vj  >  produces 

<Vp>  (Sip  +  TVi,p )  =  r{w.t;p}  <  p  >,p  <  p  >"'  . 

where  r  =  (Mtk/e)/(\  +  Mt(V/e  —  1))  The  model  is  now  a  set  of  three  coupled  linear 
algebraic  equations  of  the  form  Aip  <  Vp  >=  bi.  Inspection  of  the  equations  reveals  two 
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aigniiicaiit  features:  1)  the  mass  flux  in  one  direction,  as  might  be  expected  from  continuity 
considerations,  is  influenced  by  the  mass  flux  in  another  direction  and  2)  the  contraction  of 
the  density  gradient  on  the  Reynolds  stress  allows  coimtergradient  transfer.  In  simple  cases 
this  set  of  equations  is  easily  solved  by  hand.  Performing  the  general  inversion  the  model 
can  be  written  in  symbolic  form  as 


<  Vi  >=  TJiiivjVp]  <  p>^<  p>~^ 

where  Tij  =  {Sij  +  rK,j This  is  an  anisotropic  eddy-viscosity  model  in  which  the  eddy- 
viscosity  tensor  i/jp  =  is  a  function  of  the  Reynolds  stresses  and  the  mean  flow 

gradients.  Though  this  form  suggests  the  structure  of  the  model  it  is  not  in  a  form  most 
suitable  for  computation.  Recourse  to  the  Cayley-Hamilton  theorem  allows  Tij  to  be  written 
in  terms  of  the  invariants  and  the  first  and  second  powers  of  the  matrix: 

=A^-IaA  +  IIa1 

Substituting  A  =  1  -|-  r  W  produces  an  expression  for  the  inverse 

+  tW)-'  =  (1  -  /^  +  IIa)1  +  (2  -  U)tV\  +  T»(W)> 

and  the  final  model  can  be  written,  in  ascending  powers  of  ratios  of  time  scales,  as 

<  Vi  >=  rluoSij  -I-  i/irVi,j  Vk,j  KvjUp}  <p>,p<p  . 

The  nondimensional  "viscosity”  coefficients,  vq,  Vi,  1/2  are  known  in  terms  of  the  mean 
deformation;  they  are  not  phenomenological  parameters  that  require  calibration  to  experi¬ 
ments  which  then  limit  the  application  of  the  model  to  flows  not  too  different  from  that  for 
which  the  model  has  been  calibrated.  Only  one  phenomenological  assumption  -  to  obtain 
the  relaxation  model  for  the  correlation  with  the  fluctuating  dilatation  -  has  been  made.  The 
invariants  and  the  viscosity  coefficients  are  given  in  the  Appendix  that  summarizes  the  final 
form  of  the  model. 

5.  Discussion  and  implementation  of  the  mass  flux  model  in  simple  flows 

Formidable  as  the  algebraic  expression  for  the  mass  flux  may  appears  there  are  some 
simple  expressions  for  the  <  t;,-  >  possible.  Though  the  representation  is  valid  for  arbitrary 
three-dimensional  flows  several  cases  with  two-dimensional  mean  fields  are  investigated  in 
order  to  understand  the  effects  of  different  mean  deformations.  One  three  dimensional  field 
is  considered  in  order  to  anticipate  the  effects  the  three-dimensionality  of  the  flow  might 
have  on  the  mass  flux  expressions. 
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Case  1:  Isotropic  turbulence  with  small  velocity  gradients,  rW  <<  1 

In  this  case  rW  «  1,V  =  0  and  the  time  scale  r  =  {Mtkfe)/{1  +  Mt{Vle  —  1)).  The 

eddy-viscosity  t«isor  assumes  the  form  ~  Mt{k/e)2kSjpf{l  —  Mt)  and  the  model  is 

<  <  P>^<  P>~^  ■ 

This  can  be  compared  to  the  usual  eddy-viscosity  model:  <  Vi  >=  {ftr/  <  P>^  Pr*  p)  < 
p  >,j  in  which  pr  =  <  P>  k^/^  thus 

<  V,  >~  {k/e)k  <  p  >„•  <  p  >~*  . 

The  usual  eddy-viscosity  form  misses  the  dependence  on  Aft  which  is  necessary  if  the 
<  Vi  >  are  to  vanish  in  the  absence  of  compressibility  effects.  Thus,  apart  from  the  Aft 
scaling,  a  scalar  viscosity  assumption  will  work,  in  the  limit  of  an  isotropic  turbulence  with 
negligible  mean  velocity  gradients.  Note  that  this  form  in  a  botmdary  layer  flow  with  cross¬ 
stream  density  gradient  cannot  predict  a  streamwise  mass  flux.  It  can  only  predict  a  mass  flux 
down  the  density  gradient.  In  problems  of  engineering  interest  there  will  be  countergradient 
transport,  as  has  been  seen  in  the  Afa  =  4.5  data  of  Dinavahi  et  al.  (1993),  and  an  eddy- 
viscosity  gradient  transport  hypothesis  is  inappropriate.  These  inadequacies  have  also  been 
noted  by  Taulbee  and  VanOsdol  (1991). 

The  major  shortcoming  of  the  eddy- viscosity  assumption  is  realizability  and  its  impact  on 
computability  in  compressible  clostu^.  In  the  Reynolds  stress  equations  for  arbitrary  mean 
flow  accelerations,  a  gradient  transport  assiunption  for  <  v,-  >  can  cause  the  acceleration  pro¬ 
duction  mechanism  to  destabilize  the  computations.  For  example  if  {vaVa},  in  the  Reynolds 
stress  equations  above,  vanishes,  <  Va  >  must  also  vanish  in  order  to  keep  that  eigenvalue 
of  the  Reynolds  stress  from  going  negative,  as  a  finite  <  v,-  >  in  <  v,-  ><  p  >  D/Dt  Vj 
will  cause  negative  energies  for  arbitrary  mean  acceleration.  This  cannot  be  accomplished 
with  the  eddy- viscosity  form  of  the  model.  Gatski  (1993)  has  used  a  scalar  eddy- viscosity 
representation  and  foimd  it  to  be  computationally  destabilizing.  Zeman  and  Coleman  (1991) 
have  also  pointed  out  that  inadequate  representations  of  the  mass  flux  can  destabilize  com¬ 
putations  in  flows  when  the  acceleration  terms  is  important.  This  occurs,  for  example,  in 
the  passage  through  a  shock  or  in  flows  in  which  the  mean  strain  or  mean  dilatation  are 
important. 

Case  2:  Anisotropic  turbulence  with  small  mean  velocity  gradients,  rVV  <<  1 
In  the  case  rVV  <<  1  and  when  the  turbulence  is  anisotropic  the  expression  for  the  mass 
flux  becomes 

<  Vi  >=  T{viVp}  <p>,p<p  >-^ 

where  t  =  (Mtkfe)/{1  +  Mt{V/e  —  1)).  Note  that  this  expression  for  the  mass  flux  allows  for 
countergradient  transfer  and  is  realizable:  the  mass  flux  in  the  direction  of  the  principal  axis 
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in  which  the  eigenvalue  of  the  Reynolds  stress  vanishes  will  also  vanish.  It  is  interesting  to 
compare  the  expression  for  the  mass  flux  to  Zeman’s  (1993).  In  an  equilibrium  turbulence, 
for  which  V  =  e,  and  in  the  linait  of  small  velocity  gradients  the  present  model  simplifies,  to 
within  a  constant  of  proportionality,  to  Zeman’s  model.  Comparisons  of  this  model  with  the 
DNS  data  in  a  wall  bounded  compressible  flow  shows  that  it  does  not  successfully  capture  the 
results  known  from  DNS.  The  neglected  terms  involving  the  velocity  gradients  are  essential. 
After  all  the  terms  with  the  mean  velocity  gradients  represent  the  production  terms,  which 
are  typically  not  negligible,  in  the  mass  flux  equation.  Computations  with  the  neglected 
velocity  gradient  captiu%  the  near  wall  behavior  very  nicely,  as  will  be  seen  in  the  next  case. 
Case  3:  Simple  shear,  V,-,j  = 

In  prpblems  of  engineering  interest  the  turbulence  will  be  anisotropic  and  there  will  be 
nonnegligible  gradients  in  the  mean  velocity  field  and  the  production  terms  for  the  mass  flux 
need  to  be  included  in  the  algebraic  expression.  For  the  simple  shear,  V^,j=  V\,2  8ii8j2,  s 
surprisingly  simple  expression  for  the  mass  flux  is  possible.  The  computation  is  easily  carried 
out  .by  hand  using  <  Vp  >  {8ip  +  TVi,p)  =  r{v,t;p}  <  p  >,p<  p  >~*.  Using  the  inversion 
formula  the  invariants  of  A  are  I  a  =  3,  IIa  =  3,  IIIa  =  1  and  the  viscosity  coefficients 
take  on  the  simple  values  i/q  =  l,i/i  =  — 1,1/3  =  1  and,  as  the  square  of  the  mean  velocity 
gradient  is  zero,  the  expression  for  the  mass  flux  becomes 

<  Vi  >=  T[8ij  -  rVi,2tfrt5p]{v,i;p}  <p>,p<p  >"* 

The  streamwise  and  cross-stream  components  of  the  Favre  fluctuation  mean,  in  a  flow  with 
only  a  cross-stream  density  gradient,  become 

<  vi  >=  r[{vit;3}  -  tVu2  {V3V3}]  <  p  >,2  <P 

<  V2  >=  r{t;3t;3}  <  p>,2<P  >~*  • 

Note  that  the  effects  of  the  production  of  <  v,-  >  by  the  mean  shear,  proportional  to  V^,3, 
are  included  in  the  expressions  for  the  streamwise  mass  flux.  This  is  similar  to  the  normal 
Reynolds  stresses  in  a  unidirectional  shear:  the  production  mechanism  is  in  the  equations 
of  the  streamwise  component  of  the  energy  and  therefore  it  is  larger  than  the  spanwise  and 
cross-stream  components  of  the  turbulence  energy.  Computations,  shown  in  Figure  3,  with 
this  model  are  very  successful  for  the  streamwise  component  <  wi  >.  The  peak  in  <  ui  >  is 
captured  surprisingly  well  in  size  and  location.  This  behavior  can  not  be  captured  without  the 
inclusion  of  the  1^,3  {U3V3}  term.  The  small  velocity  gradient  limit  expression,  case  2  above, 
which  is  essentially  the  algebraic  form  both  Zeman  and  Cole  (1991)  and  Rubesin  (1990) 
substantially  underpredicts  the  near-wall  peak  of  the  mass  flux.  Note  that  a  streamwise 
mass  flux  is  engendered  by  a  cross-stream  density  gradient.  This  is  a  behavior  that  am 


isotropic  eddy-viscosity  model  cannot  predict;  such  a  model  predicts  a  zero  streamwise  mass 
flux. 

The  predictions  of  the  cross-stream  component,  <  >,  are  less  successful.  This  is  be¬ 

cause  there  are  no  large  production  terms  in  the  <  03  >  expression  to  mask  the  inaccuracies 
of  the  linear  relaxation  model  assumed  for  the  dilatational  correlation  in  a  nonequilibrium 
"newly  formed”  turbulence.  The  present  temporal  DNS  is  started  from  a  laminar  profile 
and  computed  through  transition.  The  data  shown  in  the  figures  represents  a  flow  approxi¬ 
mately  three  eddy-tumovers  past  the  transition,  £Sxf{kU)  ^  3.  The  turbulence  field  is  not 
fully  developed,  retaining  vestiges  of  the  initial  conditions;  a  lineair  relaxation  model  for  the 
correlation  with  the  fluctuating  divergence  would  not  be  expected  to  do  well  in  such,  a  more 
or  less,  transitional  flow. 

The  poor  agreement  in  the  expression  for  the  cross-stream  mass  flux  cannot  be  explained 
by  the  fact  that  the  data  comes  from  a  temporal  DNS.  The  expression  for  the  mass  flux  model 
is  from  its  evolution  equation  which  is  independent  of  the  mean  flow  equation  calculation. 
Thus,  the  pr  blem  often  seen  in  comparing  Reynolds  stress  model  calculations  to  temporal 
DNS  simulations,  in  which  there  is  a  forcing  term  in  the  equations  to  compensate  for  the 
botmdary  layer  growth,  do  not  appear  here. 

Case  4:  Plane  strain  with  mean  dilatation 

For  aplane  strain  with  arbitrary  non-zero  dilatation,  VJ,,-  =  Vi,i  SnSji+Vtti  SaSjj,  and  the  vis¬ 
cosity  coefficients  take  on  the  following  simple  values  i/q  =  1, 1^1  =  — (l+r(Vi,i  +V2,2))fIIlAi 
U2  =  If  I II A  where  IIIa  =  (1  +  tVi,i  )(1  -h  rV3,3).  The  fluxes  are  given  by  the  simple  ex¬ 
pressions 

<  <P>,p<P  >■' 

<  >=  <P>yp<P  >"*  • 

Clearly  the  model  is  fully  realizable  and  the  destabilization  of  more  rudimentary  models, 
noted  by  Zeman  and  Coleman  (1991),  for  a  flow  with  large  normal  strmn  is  not  an  issue. 
Note  that  in  very  high  strdns,  say  the  normal  passage  through  a  shock,  the  dependence  on 
the  phenomenological  parameter  r^,  absorbed  in  r  is  lost.  Here,  again,  r  =  {Mtkfe)f{l  -I- 
Af,(P/e-l)). 

In  mean  mean  field  with  a  large  dilatational  component,  it  is  more  useful  to  consider  a 
mean  velocity  gradient  described  of  the  form  Vi,j  =  Vi,i  ((1  +  —SaSji).  The  viscosity 

coefficients  take  on  the  following  values  vq  =  1,  =  —1(1  -|-  DrV^,i  )/IIIa,  1^2  =  ^IIHa- 

The  streamwise  and  cross-stream  components  of  the  Favre  fluctuation  mean,  in  a  flow  with 
arbitrary  density  gradient,  become 

<  >=  <p>,p<p  >-* 

<  '’a  >=  <  /)  >,p  <  />  >-* 
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Case  5:  Arbitrary  two-dimensional  mean  velocity  gradients 

For  an  arbitrarily  complex  two-dimensional  flow,  such  as  the  developing  wall  bounded  tur¬ 
bulent  boundary  layer  with  separation,  Vi^j  =  [Vi,i ,  Vi,2 , 0],  [Vj,! ,  Vjjj  ,  0]  is  a  suitable  rep¬ 
resentation  for  the  velocity  gradient  field.  The  viscosity  coeflScients  are  given  by  i/q  =  1, 

Ui  =  —(1  -f  TD)fIIlA,  =  1////^  where  D  =  Vi,i  -|-V2,j.  The  mass  fluxes  are  given  by 

<  ~  {^Vp}TVi,2]  <  p  >,p  <  p  >-^ 

<^>=  77^[(1  +  ){t>2Up}  -  {wiVpItVj,!  ]  <  p  >,p  <  p  >-* 

where  IIIa  =  1+tD  +  -V„2  Va,, ). 

Case  6:  Arbitrary  three-dimensional  strain  with  simple  shear 
In  a  general  three-dimensional  flow  the  expressions  for  the  invsuriants  are  somewhat  more  com¬ 
plicated.  The  simplest  case,  a  simple  shear  with  arbitrary  three  dilatation,  is  chosen.  The  ve¬ 
locity  gradients  are  represented  by  V^„  =  [V^,i ,  Vi,a , 0],  [0,  Va,a ,  0],  [0, 0, 1^,3  ].  The  square  of 
the  velocity  gradient  is  given  by  14,j  =  [(Vi,i  )*,  K,a  (Vi,i  -f  Va,2 ),  0],  [0,  (Va,a  )^,  0][0, 0,  (14,3  )*]• 
The  invariants  of  A  are  /^  =  3  -f  tD,  IIa  =  3-1-  2rD  -|-  14,a  +14,a  Hjs  +14,3 14 ,i ), 

and  I  IIa  =  (1  +  )(1  +  ^14,2  )(1  +  ^14,3).  Here,  as  usual,  D  =  14  u  is  the  mean  dilata¬ 

tion.  The  viscosity  coefficients  are  a  little  more  complicated  •  the  three-dimensionality  of  the 
flow  now  affects  the  zeroth-order  viscosity  coefficient.  In  the  two-dimensional  flows  i^o  =  1; 
here  i/q  =  (1  +  tD  +  r’d(14,i  14,a+14,a  14,3+l^r,3  Vi,i  ))////a-  The  higher  order  viscosity 
coefficients  are  given  by  i/j  =  —(1  -|-  tD)IIIa,  *^2  =  and  the  fluxes  are  written  as 

^  "i"  ^(^’2  +14,3 )  +  t*14,2  14,3  ){vit;p} 

-{v2Vp)tVu2  (1  +  rl4,3 )]  <  p  >,p  <  p 

<  wa  >=  7^1(1  +  "^(14,1  +14,3)  +  ’■*14,1 14,3){v2Vp}]  <  P  >,p<  P  >~^ 

6.  Summary  and  Conclusions 

The  fluctuating  Favre  velocity  mean,  <  v,  >,  is  the  first-order  form  of  a  second-order 
moment,  the  mass  flux,  <  pvi  >=  —  <  p  ><  t>,  >.  The  mass  fluxes  quantify  the  difference 
between  Reynolds  statistics  and  the  density-weighted  Favre  statistics,  (/,  =  K  +  <  u,  >  and 
Ui  =  Vi  —  <  Vi  >,  and  can  be  thought  of  as  measuring  the  effects  of  compressibility  due 
to  variations  in  density.  The  effects  of  the  mean  density  gradients  on  the  anisotropy  of  the 
turbulence  are  fully  parameterized  by  the  mass  flux. 

An  algebraic  representation  for  the  mass  flux  has  been  derived  from  the  transport  equa¬ 
tion  for  the  Favre  fluctuation  mean  using  the  structural  equilibrium  assumption.  The  mass 
flux  is  found  to  be  proportional  to  the  mean  density  gradients  with  an  anisotropic  eddy- 
viscosity  that  depends  on  both  the  Reynolds  stresses  and  the  mean  velocity  gradients.  The 
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appearance  of  the  mean  velocity  gradients  in  the  tensor  eddy- viscosity  reflects  their  presence 
in  the  production  terms  in  the  evolution  equation  for  <  v,-  >.  The  model  predicts  coun¬ 
tergradient  transfer  and  shows  that  mean  density  gradients  in  one  direction  can  produce  a 
mass  flux  in  a  different  direction.  It  form,  valid  for  a  general  three-dimensional  flow,  is 

<pVi>=-  <p><  Vi  >=  -T[vo6ij  -I-  I/itK„  +t/2T^9i,k  Vkyj  ]{v,Vp}  <  p  >,p 

where  t  =  {Mik/e)/{l  4-  MtiV/e  -  1)).  The  viscosity  coefficients,  are  known 

functions  of  the  mean  velocity  gradient,  given  in  terms  of  the  invariants  of  the  tensor  A  = 
1  -f  tW.  They  are  not  adjustable  "tuning”  coefficients. 

The  derivation  of  the  expression  for  <  Vi  >  has  involved  a  minimum  number  of  assump¬ 
tions  regarding  the  physics  of  compressible  turbulence.  It  is  useful,  however,  to  keep  in  mind 
some  of  the  approximations  to  2M:count  for  possible  discrepancies  and  to  anticipate  the  classes 
of  flows  in  which  the  present  form  of  the  model  may  be  inadequate.  The  assumptions  used 
are: 

1)  The  derivation  of  an  C?  (<  pp  /  <  p  >)  set  of  evolution  equations  for  the  <  Wj  > 
showed  that  the  imclosed  terms  involving  correlations  with  the  fluctuating  pressure  and 
stress  are  higher  order  effects  and  can  therefore  be  neglected.  In  the  evolution  equations 
there  is  only  one  unclosed  term,  the  fluctuating  <  ViVk^k  >  covariance. 

2)  The  form  of  turbulent  diffusion  terms  appearing  in  the  <  u,-  >  equation,  are  found 
to  scale  with  the  density  intensity,  <  pp  /  <  p  >,  for  arbitrary  inhomogeneity  and 
can  therefore  be  neglected.  This  is  consistent  with  the  truncation  of  the  equation  as  (< 
pviVp  ><  p  >“*),p=  ({w$Wp}-  <  ViVp  >),p  is  an  (<  pp  >*/*  f  <  p  >)  quantity.  The 
difference  between  {wiWp}  and  <  ViVp  >  has  been  seen  to  be  small  in  the  wall  bounded  flow 
at  Afa  =  4.5  of  Dinavahi  and  Pruett  (1993),  as  seen  in  Rjstorcelli  et  al.  (1993). 

3)  The  structural  equilibrium  assumption,  D/Dt  [<  Vi  ><  Vj  >  /{vpVp}]  =  0,  is  used 
to  produce  an  algebraic  expression  for  the  mass  flux  equation.  This  allows  the  material 
derivative  to  be  expressed  in  terms  the  production  and  dissipation  of  the  turbulence  energy. 
Pc’*  more  rapidly  varying  flows  in  which  the  structiiral  equilibrium  is  not  expected  to  yield 
results  of  adequate  accuracy  it  is  possible  to  carry  the  full  differential  equation  for  the  mass 
flux.  Near  solid  boundaries,  were  the  mass  flux  is  most  important,  a  structural  equilibrium 
is  expected  to  be  achieved  rapidly  and  the  algebrsuc  form  is  adequate.  It  is  this  fact,  coupled 
with  the  density  intensity  truncation  of  the  evolution  equation,  that  enables  the  mass  flux 
expression  to  be  used  all  the  way  to  the  wall  without  any  ad  hoc  wall  function  corrections. 

4)  The  algebraic  truncation  of  the  evolution  equation  for  <  Vi  >  involves  one  unclosed 
term,<  ViVk,k  >•  It  has  been  assumed  that  it  can  be  represented  as  a  linear  relaxation  term, 
<  ViVkjk  >=  —  <  V,  >  fti  where  =  Mtk/e.  This  model  for  the  covariance  with  the  fluc¬ 
tuating  dilatation  is  expected  to  be  adequate  for  most  quasi-equilibrium  quasi-homogeneous 
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turbulence  fields  in  whidi  the  production  terms  play  a  major  role.  For  flows  in  which  the 
pressure  dilatation  covariance  plays  a  major  role  in  transferring  energy  from  its  kinetic  to 
potential  modes  it  may  be  necessary  to  reevaluate  the  adequacy  of  the  lineu  relaxation 
model. 

The  form  of  the  mass  flmr  model  presented  does  not  include  effects  associated  with 
rotation  or  body  forces.  Both  of  these  effects  can  be  easily  incorporated  as  they  do  not 
require  any  additional  modeling;  it  is  simply  a  matter  of  retaining  the  extra  terms  in  the 
algebraic  truncation  of  the  evolution  equations.  There  is  an  exception;  at  rs|>id  rotation  rates 
the  neglected  pressure  covariance  becomes  important  and  the  truncation  of  the  evolution 
equation  used  to  obtain  the  model  is  no  longer  valid. 

The.  model  is  realizable  for  most  simple  flows  though  a  general  proof  of  its  realizability 
has  not  been  found.  These  realizability  aspects,  and  the  fact  that  the  destabilizing  pr(^>erties 
associated  with  isotropic  eddy-viscosity  models  do  not  appear  in  this  mass  flux  model,  are 
expected  to  make  it  computationally  robust. 

In  the  moment  evolution  equations  for  a  compressible  turbulence  the  mass  fluxes  appear 
in  several  places.  In  the  mean  momentum  and  energy  equations  the  mass  flux  spears  in 
five  different  locations,  Ristorcelli  (1993),  and  modeling  Ui  a  Vi  ignoring  the  contribution  of 
the  mass  flux  has  been  shown  to  be  inadequate.  In  the  Reynolds  stress  equations  the  mass 
flux  determines  the  relative  importance  of  the  production  by  the  mean  flow  accel^ation,  it 
contributes  to  the  pressure  fluxes  and  the  viscous  fluxes.  It  is  clear,  given  the  number  of 
times  it  occurs  in  the  moment  evolution  equations,  that  an  accurate  model  for  the  mass  flux 
is  necessary  for  complex  compressible  turbtdent  flows  of  aerodynamic  interest.  This  is  to 
assess  the  magnitude  of  the  mass  flux  in  various  flows  and  to  include  it  in  a  computational 
model  when  it  is  important.  There  are  classes  of  compressible  flows  in  which  the  contribution 
from  the  mass  flux  are  expected  to  be  small  and  its  inclusion  in  a  computational  model  is 
unnecessary. 

It  is  expected  that  the  mass  flux  will  not  make  much  of  a  contribution  to  usual  uni¬ 
directional  shear  flows  >.'ach  as  the  flat  plate  boundary  layer  and  diverse  free  shear  layers, 
unless  there  are  large  density  gradients.  The  mass  flux  terms  are  expected  to  be  important 
in  more  complex  flows:  these  include  flows  in  which  there  are  mean  density  gradients  due 
to  large  Mach  number  or  combustion,  separation  or  reattachment  (inflection  points),  cold 
wall  boundary  conditions,  mean  dilatation,  shocks,  adverse  pressure  gradients,  or  strong 
streamwise  accelerations  sudi  as  those  occurring  in  ramp  type  flows. 
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Appendix:  Final  form  of  the  algebraic  mam  flux  model 

The  pertineat  facts  concerning  the  computati<Hial  impionentation  of  the  mass  flux  rep¬ 
resentation  are  presented  h^.  The  gmeral  form  of  the  model,  valid  for  a  three-dimensional 
flow  in  an  inertial  system  without  body  forces,  is 

<  Vi  >=  T[v^ij  +  UiTVi,j  +MaT*V;,/k  Vk,i  <  p  >,p  <  P  >■* 

wh«re  the  viscosity  coefficients  are  functions  of  the  invariants  of,  A  =  1  -f  r  W,  and  are 
given  hy  Vo  =  {I  —  Ia  +  JIa)/IIIa,  vi  =  (2  —  I  a)!  Ill a^  and  1/3  =  If  II I  a-  The  invariants 
for  the  tensor  are  giv«t  by 

lA-i  <A> 

IIa=  1/2(<A>^-<A^>) 

.  IIIa  =  l/6(<  A>^-Z<A><A^>+2<  A®  >) 

for  which  <>  indicates  the  trace  of  the  enclosed  matrix.  The  time  scale  is  defined  as 
r  =  {Mtk/£)/(l  +  Mt{Vfe  —  1))  The  various  traces  are  straightforward  to  compute  using 
their  definition.  Their  significance  can  be  understood  when  they  are  recast  in  terms  of  the 
mean  dilatation,  rotation  and  strain: 

<  A  >ss  3  +  tD 

<A>>s  3  +  2tI?  +  t»(<5^>  +  <H^»>] 

<  A^  >*  3  +  3rX?  +  3r»(<  5*  >  -I-  <  >]  +  r®(<  > +2  <  SW^  >  +  <  >] 

D  =  is  the  mean  dilatation  and  the  strain  and  rotation  tensors  are  defined: 

«nd  Wij  s=  \l2\yi^j—Vjyi\.  The  term  tD  can  be  thought  of  as  a  ratio  of 
fluctuating  to  mean  dilatation  time  scales.  It  is  understood  to  be  order  one  or  smaller, 
tD  <\. 
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Figure  2.  The  cross-stream  second  derivative  of  the  Reynolds 
mean,  Favre  mean  and  Favre  fluctuation  mean  profiles  that  ap¬ 
pear  in  the  viscous  stress  terms  in  the  Ma  4.5  wall  bounded 
DNS  of  Dinavahi  and  Pruett  (1993). 
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